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Abstract. We investigate the observational signatures of the holographic dark energy model in this 
paper, including both the original model and a model with an interaction term between the dark en- 
ergy and dark matter. We first delineate the dynamical behavior of such models, especially whether 
they would have a "Big Rip" for different parameters, then we use several recent observational data 
to give more reliable and tighter constraints on the models. The results favor the equation of state 
of dark energy crossing —1, and the universe ends in the "Big Rip" phase. By using the Bayesian 
evidence as a model selection criterion to make the model comparison, we find that the holographic 
dark energy models are mildly favored by the observations compared with the ACDM model. 
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INTRODUCTION 

The nature of dark matter (DM) and dark energy (DE), which constitute about 95% of 
the total cosmic density, are among the most important problems in modem physics and 
astronomy [1]. Dark energy may be a problem that has to be solved in the framework of 
fundamental microscopic physics such as string theory [2]. 

Here we consider the dark matter-dark energy interaction in the case of the so-called 
holographic dark energy (HDE) model [3, 4, 5, 6, 7]. This is based on the holographic 
principle [8], which was inspired by the Bekenstein entropy bound of black holes [9]. If 
the quantum zero-point energy density Pa is relevant to an ultraviolet cut-off, the total 
energy of the whole system with size L should not exceed the mass of a black hole of 
the same size, thus we have L^pA < LM^^. The largest infrared cut-off L is chosen by 
saturating the inequality so that the HDE density is 

Pde = 3^Ml,L-^ , (1) 

where c is a numerical constant, and Mpi = 1 /\/S7rG is the reduced Planck mass. L 
should be the size of the future event horizon, since only this case would result in a dark 
energy component which drives the accelerated expansion of the universe [6], i.e. 

, , r da 

L = Reh{t)=a{t) / -r^. (2) 
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If we consider the dark matter and dark energy interaction, the energy conservation 
equation will be written as 

Pn, + 3Hp^ = Q, (3) 




Pde + 3//(l+Wde)pde = -e, (4) 

where pm, Pde and wje matter density, HDE density and HDE equation of state 
respectively. Here we consider the physically plausible interaction as 2 = 3a//pde (Of is 
a parameter describing the strength of interaction), which has been widely used for the 
interaction between the matter and massive scalar field [10]. Therefore, by combining 
the above equations and Friedmann equation, we quickly obtain the following equations 
which describe the evolution of the fractional energy density Q^ej Hubble parameter 
H{z), and also the HDE effective equation of state (EEoS) (see [7] for details): 

+ [( 1 - Qde) ( 1 + - v^) - 3 aQde] =0, (5) 
az l+z c 

^ = -:^[lt2de(l +3a + -v^) - h (6) 
az i+z 2. c I 

>Veff(z)=Wrfe + « = ---- . (7) 
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By solving the Eqs. (5), (6), and (7), we can know the dynamics of the HDE model as 
follows: For the non-interacting HDE model, i.e. a = 0, the evolution of HDE density 
and its EEoS will be affected by the value of parameter c as shown in Fig. 1. If c = 1, 
the EEoS of HDE will turn out to be — 1 (Cosmological Constant), and the universe will 
end up in a de-Sitter phase. If c > 1 , the EEoS is always greater than — 1 , and the HDE is 
just like the quintessence dark energy. On the contrary, if c < 1, the EEoS will cross —1 
at some time, and the universe will end up in a "Big Rip", indicating the HDE resembles 
the behavior of phantom dark energy [11]. 

For the interacting HDE model (IHDE), we vary the value of a but fix c = 1 to see 
the interacting effect. In Fig. 2, we see that if oc < 0, the dark matter "decays" into IHDE 
(Eqs. (3) and (4)), which will make the HDE density increase rapidly at the later era of 
cosmic evolution. Thus, effectively the IHDE will resemble the phantom dark energy 
with EEoS less than — 1 in the near future (see the left of Fig. 2). On the other hand, if 
IHDE "decays" into dark matter, i.e. a > 0, the IHDE density would decay more rapidly, 
making the IHDE behave like a quintessence field. 




FIGURE 2. EEoS and IHDE density in the interacting holographic dark energy model for different 
value of a. Left: EEoS; Right: HDE density. Here we set £2,„o — 0.27. 



However, we want to ask the inverse questions: what values of parameter c and a are 
preferred by the current observational data? In addition, which of these dynamic dark 
energy models are favored by current observational data compared with the concordance 
ACDM model? 

METHODOLOGY 

We utilize several data sets to constrain the parameters of the HDE and IHDE model, 
including the 182 high-quality type la supernovae [12], the baryon acoustic oscillation 
measurement from the Sloan Digital Sky Survey [13], 42 latest X-ray gas mass frac- 
tion data from Chandra observations [14], 27 GRB samples generated with Epeak — 
Ey correlation [15], and the CMB shift parameter from the WMAP 3 year results 
[16]. To break the degeneracy and explore the power and differences of the con- 
straints for these data sets, we use them in several combinations to perform our fitting: 
SNsei + BAG, SNsei + BAG + fgas, and SN^ei + BAG + fgas + GRB + CMB. 

For comparing different models, one must choose a statistical variable. The Xmin 
the simplest one and is widely used. However, for models with different numbers of 
parameters, the comparison using may not be fair, as one would expect that models 
with more parameters tend to have lower x^. Instead, we use the Bayesian evidence (BE) 
as a model selection criterion. The Bayesian evidence of a model M takes the form 

BE = j ^{d\0,M)p{e\M)dd, (8) 

where ^(d|0,M) is the likelihood function given the model M and parameters 0, and 
p(0 \M) is the priors of the parameters. The BE may be the best model selection criterion, 
as it is the average of likelihood of a model over its prior of the parameter space and 
automatically includes the penalties of the number of parameters and data, so it is more 
direct, reasonable and unambiguous than the Xmin model selection [17]. The logarithm 
of BE can be used as a guide for model comparison, and we choose the ACDM as the 
reference model: AlnBE = \nBE„iodei — IhBEacdm- The strength of the evidence for 
the model is considered according to the numerical value of BE: AlnBE < 1 (Weak); 



1 < AlnBE < 2.5 (Significant); 2.5 < AlnBE < 5 (Strong to very strong); AlnBE > 5 
(Decisive). We use the nested sampling technique to compute BE [18, 19]. 



RESULTS 



In Table 1, we give the best fits and the Icr CL of the HDE model parameters, as well as 
the value of In ABE for the three data set combinations. The best fit of c varies slightly 
across the different data sets, it is 0.761 for the SN+BAO data set, but decreases slightly 
when the /gas, GRB and CMB data are included. However, for all data sets, we have c < 1 
at more than 1.5a, indicating that the EEoS of HDE will cross —1 at some redshift, and 
the energy density of HDE will diverge in the future (see the blue curve in Fig. 1). The 

TABLE 1. The fitting result for the HDE model. 

I SN + BAO I SN + BAO + f.as I SN + BAO + f^as + GRB + CMB 



C 

AlnBE 



h0.020 
-0.020 
0.154 



0.273 

0.761_o.„7 
0.09 ±0.12 



-0.021 
-0.018 
0.130 



0.270 
0.63±0.18 



276+" "^^ 

^•'4°-0.093 

0.65±0.18 



HDE model fits about equally well (In BE=0.09) as the ACDM when we only use the 
SNIa and BAG data. With fgas, GRB and CMB data added, it fits mildly better than 
the ACDM, but with the data presently available the difference is not significant (In BE 
=0.63 -0.65). 

In Table 2, we show the fitting results for the IHDE model. From the best fits 

TABLE 2. The Fitting results the IHDE model. 

I SN + BAO I SN + BAO + fgas | SN + BAO + fga, + GRB + CMB 



ii;iiO 
C 

a 

AlnBE 



0.272 



h0.023 
-0.022 

-0 020+"-"*^ 
0.41+0.12 



275+''"2i 

^•^'-'-0.021 

0.667^°:?^^ 
0.068l«:?^,^ 
0.70 + 0.18 



0.281 



hO.017 
-0.017 

-0 006+'' °2i 
0.75 + 0.18 



and la CL, the fractional energy density of dark matter is still around 0.27, and the 
parameter c is always less than 1. The difference from the HDE model is that IHDE 
has the interacting parameter a here and its value is either positive or negative, but 




FIGURE 3. EEoS and energy density evolution for the best fit IHDE models. 
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FIGURE 4. The contours of c vs. a in the IHDE model. 



always very small and around zero. If you look closer at it, zero is always covered 
in the la CL, indicating that there is little evidence for the interaction. Furthermore, 
the Bayesian Evidence is always positive but less than 1, if you compare them with 
the previous Bayesian Evidence of non-interacting case, you can see that the number is 
slightly greater than the previous case, which means that the IHDE model is a little more 
favored than the HDE model but the evidence is still not strong. 

Using the best fitting values for the three data sets, we also plot the effective equation 
of state and dark energy density in Fig. 3. We can see that all three data sets suggest the 
equation of state crosses —1 and the universe ends in the "Big Rip" in the future, just 
like the HDE model. 

The contour map for c and a is plotted in Fig. 4. In Fig. 4, the contours corresponding 
to the data sets SN+BA0+/^a.5+GRB+CMB are much tighter than the other data sets. 
This is because we use the data set GRB and CMB, which has very large redshift 
distribution so they break the degeneracy of the parameter c and a. We also mark the 
a — c values for which Weff = — 1 as a dashed line. This forms the dividing line between 
quintessence-like and phantom-like behavior. For the best fit parameters of all three data 
set combinations, c ~ 0.6, so c < 1. The value of a varies more, but all consistent with 
being within 1.5c7, which suggests the evidence for the interaction is very weak. In any 
case, for all three data set combinations, the best fit value resides in the phantom-like 
region, although a large area of quintessence-like region is also allowed. 

CONCLUSION 

In this paper we first introduced the holographic dark energy model with the interact- 
ing term Q = 3a//pde, so the non-interacting case could be viewed as the special case 
with a = 0. We illustrated the dynamical behavior of these models by choosing some 
representative values of the parameters c and a. The condition for the model to have 
a "Big Rip" is determined. Second, we utilize several data sets from the recent obser- 
vations to constrain the models. The best-fits for the three data sets are given in Table 



1 for the HDE model, and Table 2 for the IHDE model. For both the HDE and IHDE 
models, the data favors "phantom" behavior slightly, i.e. the dark energy initially has 
Wgff > — 1, but eventually crossing the phantom dividing line, and the model ends with 
a "Big Rip". However, quintessence-like behavior is also allowed with the present data. 
Next, we utilize the Bayesian evidence (BE) as a model selection criterion to compare 
the holographic models with the ACDM model for the three data sets. Both the HDE 
and the IHDE model are mildly favored by the current observational data set, although 
the evidence is weak. 

In brief, we conclude that according to the observational data, the holographic dark 
energy model, especially the interacting holographic dark energy model is mildly fa- 
vored by the current data, and for the best fit model the equation of state for both the 
HDE and IHDE crosses —1, for which the Universe ends up in a "Big Rip". 
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